{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "读取数据"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "a=79485\n",
    "b=106339"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pylab as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "(time_stamp,indicated_airspeed,true_airspeed) = np.loadtxt(open(\"d:\\\\flightdata\\\\newfile\\\\log_301_2020-3-12-14-34-46_airspeed_0.csv\",\"rb\"),skiprows=1,delimiter=\",\",usecols=(0, 1, 2),unpack=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "ref_time = (time_stamp - time_stamp[a])/1000000"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deXgUVdb48e/JRtiR1aCsbsgiESK7Aq/ooCiKM4774DKgDqKIryOKCvq6zs9tGHDBBXRExhlGBZRxBAEXFJFIkLATBAlEliCEkK07fX5/VHeTQBI6CUl3JefzPP10pbqW013p07du3bpXVBVjjDG1R1S4AzDGGFO9LPEbY0wtY4nfGGNqGUv8xhhTy1jiN8aYWsYSvzHG1DKW+E2tJCI3iMhn4Y7jeEREReT0cMdhahaxdvzGRC4RUeAMVd0S7lhMzWElfmOOIiIx4Y7BmKpkid/UaCIyQUTSROSQiKwTkRH++TeLyNdFllMRGSMim4HN4nhRRPaIyEER+VFEuvqXnSkir4rIQv92vxCRdkW21cn/2n4R2Sgivy/yWh0ReU5EfhaR3f7t1C3y+v0ikiEiu0Tk1mr5kEytY4nf1HRpwPlAY+Ax4F0RSShl2SuB3kBn4GLgAuBMoAlwDZBZZNkbgP8DmgMpwCwAEakPLATeA1oC1wEvi0gX/3rP+reZCJwOnAI86l93KPC/wEXAGcCQSr1zY0phid/UaKr6L1Xdpao+VX0f2Az0KmXxp1V1v6rmAh6gIdAJ51rYelXNKLLsJ6r6parmAxOBviLSBrgM2KaqM1TVq6o/AP8GficiAowC7vXv5xDwFHCtf5u/B2aoaqqqHgYmn8jPwpgAq8s0NZqI/AEYD7T3z2qAU0ovLGHxHYEJVV0sIlOBaUBbEfkQ+F9VzSph2WwR2Q+0BtoBvUXkQJHtxgB/B1oA9YBk5zfACRGI9k+3BpKLrLe9XG/WmBBZid/UWP5699eBu4BmqtoESMVJtiUp1sRNVaeoak+gC071zP1FXm5TZD8NgKbALpwfhC9UtUmRRwNVvRPYB+QCXYq81lhVG/g3lVF0u0Dbir1zY8pmid/UZPVxkvleABG5Begayooicp6I9BaRWOAwkEfxs4RLRWSAiMTh1PV/p6o7gI+BM0XkJhGJ9T/OE5GzVdWH80P0ooi09O/nFBH5jX+b/wRuFpHOIlIPmFTZD8CYkljiNzWWqq4Dnge+BXYD3YBlIa7eCCdJ/4pT5ZIJPFfk9fdwEvN+oCfOxV789fYX49Tb7wJ+wbmgW8e/3gPAFmC5iGQBi4Cz/Ov+B3gJWOxfZnE537IxIbEbuIwpJxGZCaSr6sPhjsWYirASvzHG1DKW+I0xppaxqh5jjKllrMRvjDG1jCtu4GrevLm2b98+3GEYY4yrJCcn71PVFkfPd0Xib9++PStXrgx3GMYY4yoiUuLd31bVY4wxtYwlfmOMqWUs8RtjTC3jijr+kng8HtLT08nLywt3KCaM4uPjOfXUU4mNjQ13KMa4hmsTf3p6Og0bNqR9+/YU6eLW1CKqSmZmJunp6XTo0CHc4RjjGq6t6snLy6NZs2aW9GsxEaFZs2Z21mdMObk28QOW9I39DxhTAa5O/MYYE+l8viPThSWN+xYGlviNMaYKNWsGa9c60zExsGNH2ctXB0v8ldCvX79yLb906VIuu+wyAObNm8czzzxTof0+9dRT5V5n5syZ3HXXXeVerzJxnkjbtm2ja9eQBs8yYXL4cOSUaCPJgQOwcSME+sNMTQ1vPGCJv1K++eabCq87fPhwJkyYUKF1K5L4K6q0OL1eb7XFYNxh8GB46aVwRxFZ7rnHed64EQYMcKY3bAhfPAGubc5Z1Lhx40hJSTmh20xMTOSl4/wXN2jQgOzsbJYuXcrkyZNp3rw5qamp9OzZk3fffRcR4dNPP2XcuHE0b96cHj16BNedOXMmK1euZOrUqezevZs77riDrVu3AvDKK6/Qr18/rrzySnbs2EFeXh733HMPo0ePZsKECeTm5pKYmEiXLl2YNWsW7777LlOmTKGgoIDevXvz8ssvEx0dzYwZM3j66adJSEjgzDPPpE6dOqW9FebPn88TTzxBQUEBzZo1Y9asWbRq1apYnDfffDNNmzZl1apV9OjRg+HDh3OP/z9bRPjyyy9JTk7m0UcfpVmzZmzcuJELLriAl19+maioKD777DMmTZpEfn4+p512GjNmzKBBgwYkJyczfvx4srOzad68OTNnziQhIYHk5GRuvfVW6tWrx4DAt8ZErF27YP/+cEcRWVasgN274YornM9n9WqYOjXcUVmJ/4RZtWoVL730EuvWrWPr1q0sW7aMvLw8Ro0axfz58/nqq6/45ZdfSlz37rvvZuDAgaxevZoffviBLl26APDWW2+RnJzMypUrmTJlCpmZmTzzzDPUrVuXlJQUZs2axfr163n//fdZtmwZKSkpREdHM2vWLDIyMpg0aRLLli1j4cKFrFu3rsz4BwwYwPLly1m1ahXXXnstf/nLX0pcbtOmTSxatIjnn3+e5557jmnTppGSksJXX31F3bp1AVixYgXPP/88a9asIS0tjQ8++IB9+/bxxBNPsGjRIn744QeSkpJ44YUX8Hg8jB07ljlz5gQT/cSJEwG45ZZbmDJlCt9++21FD4upRvHxkJsb7igiS7Nm0KIFLF8Ojz0G3brBzz+HO6oaUuI/Xsm8OvTq1YtTTz0VcM4Wtm3bRoMGDejQoQNnnHEGADfeeCPTp08/Zt3FixfzzjvvABAdHU3jxo0BmDJlCh9++CEAO3bsYPPmzTRr1qzYup9//jnJycmcd955AOTm5tKyZUu+++47Bg0aRIsWTo+s11xzDZs2bSo1/vT0dK655hoyMjIoKCgo9Yaoq6++mujoaAD69+/P+PHjueGGG7jqqquC779Xr1507NgRgOuuu46vv/6a+Ph41q1bR//+/QEoKCigb9++bNy4kdTUVC666CIACgsLSUhI4ODBgxw4cICBAwcCcNNNN/Gf//yn1PhN+LVuDRkZ4Y4i8gRaHHfu7Ew3auScGTVtGr6YakTijwRFq1Gio6ODdeAVbWe+dOlSFi1axLfffku9evUYNGhQiTcqqSojR47k6aefLjb/o48+Kte+x44dy/jx4xk+fHiw6qok9evXD05PmDCBYcOGsWDBAvr06cOiRYuAY9+ziKCqXHTRRcyePbvYa2vWrKFLly7HlOoPHDhgbfRdpnFjOHQo3FFEpl27oFUrZ7pvX0hOBn9ZJyysqqcKderUiZ9++om0tDSAY5JewIUXXsgrr7wCOCXerKwsDh48yEknnUS9evXYsGEDy5cvDy4fGxuLx+MJrjtnzhz27NkDwP79+9m+fTu9e/dm6dKlZGZm4vF4+Ne//lVmrAcPHuSUU04B4O233w7p/aWlpdGtWzceeOABkpKS2OC/arVixQp++uknfD4f77//PgMGDKBPnz4sW7aMLVu2AJCTk8OmTZs466yz2Lt3bzDxezwe1q5dS5MmTWjcuDFff/01ALNmzQopJmPC6fPPIfBV8/mOlPYTEiDKn20j4czIEn8Vio+PZ/r06QwbNowBAwbQrl27Epf761//ypIlS+jWrRs9e/Zk7dq1DB06FK/XyznnnMMjjzxCnz59gsuPHj2ac845hxtuuIHOnTvzxBNPcPHFF3POOedw0UUXkZGRQUJCApMnT6Zv374MGTKk2IXlkkyePJmrr76a888/n+bNm4f0/l566SW6du1K9+7dqVu3LpdccgkAffv2ZcKECXTt2pUOHTowYsQIWrRowcyZM7nuuus455xz6NOnDxs2bCAuLo45c+bwwAMP0L17dxITE4OtpWbMmMGYMWPo27dv8PqBiWwiR5p0+msva43Dh+HFFyFw8p2bC/XqKU899RRr164lNzeXUaNGcfDgRkq53Fd9VLVKHkAbYAmwHlgL3OOfPxnYCaT4H5ceb1s9e/bUo61bt+6YeSb8lixZosOGDavWfdr/QmQ4+eTvtVu3LzUtTTUrSxVUDx0qedncXNXzzjvy95Il1RJilZo61XnPXbqobtyoumOH6nXXHVZAW7durUuWLFFAe/W6UceNq56YgJVaQk6tyhK/F7hPVc8G+gBjRKSz/7UXVTXR/1hQhTEYY6rJL79ksGbNPCZNgkAjso0bnef8/OI3d+3cCd9/70z7fM49AG538KDz/MsvcNZZ8M9/wnffLQZg165dDPa/ye3bl9fcqh5VzVDVH/zTh3BK/qdU1f5MaJ588kkSExOLPZ588skTtv1Bgwbx8ccfn7DtGXc4cOCAf2o66enwyitw++1Hblbq2xc++ujI8jt3Os9e75EuDAJ3tpbk8OGyX68u6enOnbgAR9/Qvm2b83zWWc7zN9/ksnXr6mO2cfBgOtnZ4X0z1VLHLyLtgXOB7/yz7hKRH0XkLRE5qZR1RovIShFZuXfv3uoIs1aYOHEiKSkpxR6BdvPGVFSgcQFkUbeuc9PS7bfDjTfCp5863RSkpx9Zftcup93/nj1Hzgpyckrffrt28N//Vln4IfvNb+DRR53pBx8s3oopI8M50/n6a/jwQ0hOViCbCy+8MLjM2LFjycvLIycnvDc8VHniF5EGwL+BcaqaBbwCnAYkAhnA8yWtp6rTVTVJVZMCbdGNMZEpOzsXcEqxe/b4aNkS/PchMn26k7gDJWJwEv9JJ8GCBbB5M7Rv7/wIgFOynzYNirZezsyEtDSnGeQXXzjzsrKcO2Gr0759TtyBHjeLlklV4eyznQvcl10GO3fWAbK5/PLLAYiJieHGG28E4Ndfw3uLc5UmfhGJxUn6s1T1AwBV3a2qharqA14HelVlDMaYqnfoUAGQD0BmppdGjSAuzqnXLyx07mB97z2nu4ING+C++5wS85o1Tkm5Vy+4916nb5tnnoG77nJeCzjlFOfHYs4c+O1vnXlffgmJic6ymZmwbNmxcbVocaSKqCJVRevWOdVUATExzvWKQNcUgbOUggIoOvpnTAx4PNFANhdccAEAI0eOJCkpiZYtW3LgwC4OHy5/PCdKlSV+ce6+eRNYr6ovFJmfUGSxEUAE9FVnjKmMQ4fyCST+8eNTuPXW/Tz99NPk5eXQubPTTUHHjjB2rFMqBvj9751S/t69cO65MHcupKTAxx/Ds8861wFatXJ+FDp2dBK/zweBewgDpe1p02DECLjhhuIx7d/vlND37HHa1/fsWfZ7UIXs7OLdTnz9tdPdQqtWzhlI27bOj8zu3c7rgcS/Y4fz2rx583jvvfeKbDWbhIQEfv75Z1599VWioqI47bTT8Hp3hPUCb1WW+PsDNwH/IyIp/selwF9EZI2I/AgMBu6twhiMMdXg8OECwLmpMDv7c5YsmclDDz3E66+/zm23wauvwrffOs8Ao0ZlcMcdI9i/38vevUrjxk5F/7ffOlU4XbvCli1O0v7+e+fsYd8+J/l36eL8KGRkOKX+uXOdljSBqqVTT3U6R5s3z/nB+Ogj5yxj1Srn9bffLl76f/dd5zkqCho2hIcfdv7+7rsjZx179sD27c6Px759TuJv3ZpgqX37dmjSJIsrrriCG4r9AmXTpEkT2rRpQ0yM01FC27Ztyc1NY9euE3gAyqkqW/V8raqiqucUbbqpqjepajf//OGq6srePQ4cOMDLL78cln1bH/km0uTkeIACAB566KFgv1Bbtmzh9NOhf/9McnJyuP12Z/k9ex7no48+4sCBA6Sl7eBPfzqbmTOXct55MHu2U7UzZw783/85F0rvvttJ1r/+Ci1bOsk9I8N5Hj4cAt1Qeb3Oj8LmzfD++/DQQ3DHHTBjhlPdtG8f3HwzrFzpXHSePx9uuglOO+3Y99Snj9M6adYsp7npli3Och6Pk/g7dnTOEJz3Ce+/f6S79BUrVgAQFZV7TK+4bdu25eDB9ezcGb6WPXbnbgWVlvgLq2EkCusj30Sa7GynxB+4w/q1114DYP369RQWFtK8eXNGjBgBOPXhublOF+Reby579x4EFI9nC8uWOaX9U091StwjR8KVVzrJPeDFF2HgQKf0n1Ck4jguzqmT79nTaU3UrBl07+68tmqVc8H1+uuhXz/45ht4/HFnu506wdatsHChU7LfscOp1qlXz7k+cf31ztlEcjKcfLKzvd274fTTnbb6vXo5ZwYbN84JxtK7d2/ga+rWzTmmz6mTTz4Zr3cbmzYd2/dWdbHEX0ETJkwgLS2NxMREzjvvPAYPHsz1119Pt27djikJP/fcc8FOz9LS0hg6dCg9e/bk/PPPD/ZvU5L58+fTu3dvzj33XIYMGcJuf8Vi0dG0br75ZsaPH8/gwYN54IEH+OKLL4Lt888991wOHTrE0qVLueCCCxgxYgSdO3fmjjvuwOdvlvDZZ5/Rt29fevTowdVXX022vwiTnJzMwIED6dmzJ7/5zW/I8FdIJicn0717d/r27cu0adNO+Odq3ClQ1fPss88Wm//555+TmJgIOP9r4FwEjfJ3XBMXt5dDhxoBTo+tAc2aOe3l27Rx/h48eDDbtm0gL89pVdO1q5Pk4+OP7Ou005xEf9tt8Mc/OlU4Z57pXFPw+WDMGCe5z57tlPQDvYB8+qmT7IcMcS4G/+tfkJQEl17q/FiAcwYybx5s2zafpUs9vPZaIV9/PZN333XOHlJTC4CtR30q59OoUTpHO/nkk4G1rFpVcMxr1aXG9M45ejQntM6sdWunGVppnnnmGVJTU0lJSWHp0qUMGzaM1NRUOnTowLai7daOiXM0r776KmeccQbfffcdf/rTn1i8eHGJywb6yBcR3njjDf7yl7/w/PPHtn4N9JEfHR3N5ZdfzrRp0+jfvz/Z2dnE+78ZK1asYN26dbRr146hQ4fywQcfMGjQoGAf+fXr1+fZZ5/lhRde4MEHH2Ts2LHMnTuXFi1a8P777zNx4kTeeustbrnlFv72t78xcOBA7r///nJ9pqbmys31AgUkJiYF502cOJEnn3yS1CJjDaoqIsJ+f7OY5OT5+HzjANi3bx/gnDW/+eabDB06lMaN2wYHO4IFHDzYgW+//YHExL5Mm+Z0Jx7oDvzhh53kPWSI06oGnC6QP/nEqZZJSnJK9m3bOhd7ExKcsXCP7kLr7rthyhQYN875AQEn8ScnQ3LycOALNmy4AHBuTx42DHbvPtKctV27dmzfvh0o3pttQKtWrYC9/PJL+M7Qa0ziLytJV4devXqV2od9QHZ2Nt988w1XX311cF5+fn6py1sf+cYtcnK8gIdOnToF5/35z39m+fLlfP7558F5AwcO5H/+53+CdeAeTz5RUU7CnDRpEtu2bSM1NZXvv/+e3/72t8yZMyfYuy3sID1d6NevH5995gWiadOmDc8++yyjRo0iMzOToUNPD+7r8OHDPPLII9x//59ZvvxkRCDwFfJ4nPb2/q9NMX/9K0yc6HQz/euvv7Jhwwbq1+9LnTqK83UdiJPkX2fQoALi41+icWPnh2zVqlUkJiYyceJEnnrqqWDBq6iT/fVFBQX5qB7pwbM61ZjEH25Ff9ljYmKCVSlAsB99n89HkyZNQh4m0vrIN26Rl1cIeKhfvz6XX3457du3p1GjRixatKjY/8xXX33FV199VWTNKHy+IyXfGTNmBKf3799Pq1atitwV/AawEIDt21MBpwL/gQce4IEHHgCc71hgf7Nnz+bFF1+kUaNG9OrVi8mTVwS/QzHHyXwtWzrPt9wyhtmzZ/PWW+/w/vtncOWVzvwRI5bw4YcHOPXUTKZMgddf/w+LFhEsaLXx11EdKmGAgkDij4rKIiurNf5xl6qV1fFXUMOGDUs8qEDwnzUzM5P8/Pxg3zWNGjWiQ4cOwb7xVZXVZdx6aH3kG7fIz/cBHmJjY5k3bx5TpkwJvnb77bfTvn17XnjhhWLrOHexzgZuLXGbS5YsKZL0AbJxOvqF3NyviI1NOGadPXv2MGjQID788MNg996bNm1i2LBhPPbYYxw+fJiffvqp2DoFBQUUFBQwfvx4rr/++mIj+gUKRbfe+gc++ui14Pz27ecDzqBLJ50Ee/f+RJ06dYIj5AV+AH799ddjYmzatClxcXFERe1m6dIS37o/riN3CJ9olvgrqFmzZvTv35+uXbseU9cdGxvLo48+Su/evbnsssuKnf7OmjWLN998k+7du9OlSxfmzp1b6j6sj3zjFl99dR5QEGyrXtTUqVPZuHEj9957b7Aw0a1bN2666SacC6LzufTSS2ndunVwnZEjR5a5v7vvHovH8wvPPfccjz/+eHD+rFmz+OKLL7jqqqv4+9//DhQfAGnSpEl07NiR1157jQ0bNvDaa69Rp04dGjVqxIsvvsjs2bO59957Szy7DXwHANaudX6ADvsb8geuNQTWCRTYsrKyjoldRGjdujWFhTu58krnhrCSjBwJVTaqbEl9NUfaw/rjr5xw9JFfnex/IfxANTp6/HGX8/l8+uijj2paWpqmpqYqTmW5zp8/X1VV33nnHd2zZ48++OCDCuhVV12lq1ev1jfffDO4bNHHzJkzVVX1xx9/9Pd136vE5Y5+DBkypMzXFy5cWObrrVu3VkCHDh2qmzdv1n79+unAgQOD7/OXX34JLluSfv36aWLivQqqKSlH5ufkqHq9qq+8otq7t2oJqa9cKKU//rAn9VAelvgrxxK/qUo+n5NJ6tQp3+giBw4cCCbHFStWFHtt+fLlOnz4cN27d6+qqh48eDC47Pfffx+c/uSTT1RVNScnp8QEPWzYsJB+CAAdM2aM3nbbbQpoTExMcP6aNWu0Tp06x11/7NixwfgLCwsV0FGjRpX43q+99lpt3769/u1vqgsWOPMyM53P8bPPnOff/U5169ZyfaTHKC3xW1VPBLA+8o2bBZrfR0eX707URo0aBadbBq6m+vXu3Zu5c+cGqziLLhu4cAoEX69bt26wF8yLL76Yl156CZ/PV+y+AuemKvjDH/5QbF/jx4/n+++/Z+rUqfzud78DjtwMedZZZ9GpU6dgtebwoneSHWXAgAHB6aioKNasWVNi82uAc889l23bttGgQRa7djn9Dt13n3OPwsUXO10/T516pBXSCVfSr0GkPUor8ft8vsr9HBrX8/l8VuIPs19/dUqoDRsev6rnaPhLyzk5OSEv6/F4gtNpaWnB1wsKCnTnzp3HrPf0008roOvXr9dt27ZpWlpasSodr9cbXDYnJ0c7duyogPbr1++YfZdVBXTw4MGQ3/fixYsV0L/+dbk6nVE4j7lzVTdtCnkzx0VNK/HHx8eTmZnp1FeZWklVyczMLLGttKk+gX7zo6PrVXgboTQUePDBBzn99NOLXUAueqYQGxtb7AJxwH333cfatWvp1KkT7dq1o2PHjsEEOGTIkOA9MIE40tLS2Lp1a7F7VAJnE4MHD+att94CnFL7G2+8AUCLFi2KnZUcz+mnO/cbREX9GJyXnQ2XXw5nnBHyZirMte34Tz31VNLT07HRuWq3+Pj4YNM5Ex65uQpIhRL/I488QnJyckjLPvXUUzz1lNMR2scff8w//vGPEu+MPVpsbCydO3c+7nJFHX2z5BtvvMGuXbuIjo7m2muvJSUlhYcffpgWLVrQo0ePYBcUoWrdujUxMTHs2vUTPp/TPUUIb+WEETeUmJOSknTlypXhDsMYU4JJk2bz+OOJwCWobgt3OK7RsGFDsrOzWbx4cXAg9qImT55M69atGT16dIX3ISLJqpp09HzXVvUYYyLDjz9uAuYA28MdiqsEOkQMVB0F+Hw+fv75Zx577DFuv/12cnNP/Pi8lviNMZWiWgcIXxfDbhXo1j3WP2bjkiVL6N27NxMmTKBdkZ7jqqI/LNfW8RtjIoNqPGDX2srrzjvvZObMmaSnp5OVlcWYMWNYv359sAM7gDlz5jBs2LATvm9L/MaYSrESf8W1adOGf//73zQ+qqe2k08+mYULF1bZCHdW1WOMqRSfLw5L/BXj9M1/xIQJE9i4cSMZGRlVOqyplfiNMZXiNOfMpUmTJuEOxXXuvPPOYF1/RkZGsMvmqmaJ3xhTKbm5PiCPdevWhTsU1+natWvwZrLqHOvCqnqMMZWSmwtt2rQkIeHY/vFNaKp7gCNL/MaYSsnLE6zXDHexxG+MqZT8/ChsTB53scRvjKmU/Pxo6taN/K5fzBGW+I0xlVJQEE29etVbR20qxxK/MaZSPJ4Y6te3VOImdrSMMZXi9Vrid5sqO1oi0kZElojIehFZKyL3+Oc3FZGFIrLZ/3xSVcVgjKl6Xm8cDRpY4neTqjxaXuA+VT0b6AOMEZHOwATgc1U9A/jc/7cxxqW83lgaNowNdximHKos8atqhqr+4J8+BKwHTgGuAN72L/Y2cGVVxWCMqXqFhT7q1bOG/G5SLednItIeOBf4Dmilqhng/DgALUtZZ7SIrBSRlTa8ojGRyePxABrSmLkmclR54heRBsC/gXGqmhXqeqo6XVWTVDWpRYsWVRegMabC8vLyALEB712mShO/iMTiJP1ZqvqBf/ZuEUnwv54A7KnKGIwxVScwLKCV+N2lKlv1CPAmsF5VXyjy0jxgpH96JDC3qmIwxlQtJ/FbVY/bVGW3zP2Bm4A1IpLin/cQ8AzwTxG5DfgZuLoKYzDGVCEr8btTlSV+Vf0aKO0+7gurar/GmOpjid+d7K4LY0yF5eTkAmKJ32Us8RtjKuzQoTygwBK/y1jiN8ZUWFaWB8i1xO8ylviNMRV28GABkGvt+F3GEr8xpsIOHfJiJX73scRvjKmwrCxL/G5kid8YU2FOiT/HEr/LWOI3xlRYdnYhVuJ3H0v8xpgKy8lR7OKu+1jiN8ZU2OHDPmJivDhdcxm3sMRvjKmwnBwlLs4b7jBMOVniN8ZUWE4OxMUVhjsMU06W+I0xFZabC/HxGu4wTDlZ4jfGVFhenljidyFL/MaYCsvLi7LE70KW+I0xFVZQEE29euGOwpSXJX5jTIU5id+acrrNcUfgEpGWOMMotgZygVRgpar6qjg2Y0yEs8TvTqUmfhEZDEwAmgKrgD1APHAlcJqIzAGeV9Ws6gjUGBN5PJ44GjSwigO3KavEfykwSlV/PvoFEYkBLgMuAv5dRbEZYzQgGDoAABFfSURBVCKcx1OXRo3s4q7blJr4VfX+Ml7zAh9VSUTGGFfIzgavV2jUqEG4QzHldNxzNBG5R0QaieNNEflBRC6ujuCMMZFr6FDF6x1CgwaW+N0mlMq5W/31+BcDLYBbgGeqNCpjTMQ7fNip4qlfv36YIzHlFUriD1yyvxSYoaqri8wzxtRSMTFOwz4r8btPKIk/WUQ+w0n8/xWRhoA15TSmlsvPdzpns8TvPqUmfn/LHYDbcJp1nqeqOUAcTnWPMaYWS0v7AOhqid+FymrOuVxE0oFPgU9V9QCAqmYCmdURnDEmcuXkFALpNG/ePNyhmHIqqzlnkoi0Ay4BXhKRU4Cvgf8AX6hqfjXFaIyJSPFAHq1atQp3IKacyqzjV9Xtqvqqql4J9APmA0OAr0Tkk7LWFZG3RGSPiKQWmTdZRHaKSIr/cemJeBPGmHCoA+Rb4neh4/bVE6CqHmCx/4H/DKAsM4GpwDtHzX9RVZ8rR4zGmIjkNO5r2rRpmOMw5RXKDVyXicgqEflVRLJE5JCIZKnqzrLWU9Uvgf0nLFJjTMTwep1xdh966CGioqyvHrcJ5Yi9BIwEmqpqI1VtqKqNKrHPu0TkR39V0EmlLSQio0VkpYis3Lt3byV2Z4w50fLy8gBo0qRJmCMxFRFK4t8BpKrqieiJ6RXgNCARyACeL21BVZ2uqkmqmtSiRYsTsGtjzIkSSPzx8fFhjsRURCh1/H8GFojIF0CwJY+qvlDenanq7sC0iLwOfFzebRhjws8Sv7uFkvifBLJx2m7FVWZnIpKgqhn+P0fgDOpijHGZ/Px8QCzxu1Qoib+pqpa7N04RmQ0MApr7bwSbBAwSkURAgW3A7eXdrjEm/JwSv1KnTp1wh2IqIJTEv0hELlbVz8qzYVW9roTZb5ZnG8aYyGRVPe4WysXdMcCnIpJbtDlnVQdmjIlclvjd7bglflVtWB2BGGPcIycnD/BZ4nepsnrnbF/Wiv4RuU490QEZYyLf4cMeoMASv0uVVeL/fyISBcwFkoG9OC17TgcGAxfiXLBNr+ogjTGR5dChAiDfLu66VFm9c14tIp2BG4BbgQQgB1gPLACeVNW8aonSGBNRrMTvbmXW8avqOmBiNcVijHGJ7GwPkG+J36WsdyVjTLnl5HixxO9elviNMeXmjL5lVT1uZYnfGFNugRK/Xdx1p1Lr+EWkR1krquoPJz4cY4wbBEr8lvjdqayLu4Euk+OBJGA1zpA75wDfAQOqNjRjTKTKySkkOtqLiIQ7FFMBpVb1qOpgVR0MbAd6+PvG7wmcC2yprgCNMZEnN1eJjfWFOwxTQaHU8XdS1TWBP1Q1FWcgFWNMLZWXZ4nfzULpnXO9iLwBvIvTnfKNODdxGWNqKUv87hZK4r8FuBO4x//3lzhDKBpjaqncXCUuzhK/W4XSO2eeiLwKLFDVjdUQkzEmwuXnQ1zciRiG24TDcev4RWQ4kAJ86v87UUTmVXVgxpjIVVAgWEtO9wrl4u4koBdwAEBVU4D2VRiTMSbCFRREWVWPi4WS+L2qerDKIzHGuIbHE2UlfhcL5eJuqohcD0SLyBnA3cA3VRuWMSaSeTxCXFy4ozAVFUqJfyzQBcgH3gMOAuOqMihjTGTzeoX4eOvqy61CadWTA0wUkadU9XA1xGSMiXBer1X1uFkorXr6icg6/DdtiUh3EXm5yiMzxkQsJ/FHhzsMU0GhnKu9CPwGyARQ1dXABVUZlDEmshUWRllVj4uFdORUdcdRswqrIBZjjEt4vVHEx1uJ361CadWzQ0T6ASoicTiteqyvHmNqscLCaOrWtcTvVqGU+O8AxgCnADtxeuYcU5VBGWMil6ri9UbRoIFd3XWrUFr17ANuqIZYjDEukJOTA9SneXNL/G4VSquejiIyX0T2isgeEZkrIh1DWO8t//KpReY1FZGFIrLZ/3xSZd+AMaZ6HTp0CBAaN24Y7lBMBYVS1fMe8E8gAWgN/AuYHcJ6M4GhR82bAHyuqmcAn/v/Nsa4iJP4oVGjRmGOxFRUKIlfVPXvqur1PwIDspRJVb8E9h81+wrgbf/028CV5YrWGBN2WVlZADRsaCV+twqlVc8SEZkA/AMn4V8DfCIiTQFU9ejkXpZWqprhXy9DRFqWN2BjTHilpfkAoV69uuEOxVRQKIn/Gv/z7UfNvxXnh+C49f0VISKjgdEAbdu2rYpdGGPKSRWuueY84GNiY2PDHY6poFBa9XQ4gfvbLSIJ/tJ+ArCnjP1OB6YDJCUl2VA/xkSAwiK3blrid69QWvVcLSIN/dMPi8gHInJuBfc3Dxjpnx4JzK3gdowxYfDDD2v8U0JMTCgVBiYShXJx9xFVPSQiA3D67HkbePV4K4nIbOBb4CwRSReR24BngItEZDNwkf9vY4xLPPbYk/4ptRK/i4Xykx04uRsGvKKqc0Vk8vFWUtXrSnnpwhBjM8ZEGK9XgtOW+N0rlBL/ThF5Dfg9sEBE6oS4njGmxgkke48lfhcLJYH/HvgvMFRVDwBNgfurNCpjTIQKJPt8S/wuFuoIXB8U+TsDyKjKoIwxkUk1kDLy7eKui1mVjTEmZEUTv5X43csSvzEmZCKBZO+zxO9ilviNMeUQ538WS/wuZonfGBOyqKgjffBbHb97WeI3xpRDoJRvJX43s8RvjAlZVJRV9dQElviNMSHz+QLVO1GISJnLmshlid8YEzKPJ9BRriV9N7PEb4wJ2ZG+eizxu5klfmNMyJwS/wbgy3CHYirBEr8xJmQeD8CTwMzwBmIqxRK/MSZkHo8AHmvD73KW+I0xIfN6FfCwbNmycIdiKsESvzEmZB4PnH32WfTq1SvcoZhKsMRvjAlZYSHExlqLHrezxG+MCZnXC9HR4Y7CVJYlfmNMyAoLxUr8NYAlfmNMyHw+sAY97meJ3xgTssJCISbGSvxuZ4nfGBMyq+qpGSzxG2NCVlgoREfr8Rc0Ec0SvzEmZD5fFLGxljbczo6gMSZkPp/Yxd0awBK/MSZkPp/V8dcElviNMSGzqp6awY6gMSZkVtVTM4TlEIrINuAQUAh4VTUpHHEYY8rHSvw1Qzh/uwer6r4w7t8YU04+X7TV8dcA9tNtjAmZzyfExVnacLtwHUEFPhORZBEZXdICIjJaRFaKyMq9e/dWc3jGmJKoRltVTw0QriPYX1V7AJcAY0TkgqMXUNXpqpqkqkktWrSo/giNMcX4fD4gyqp6aoCwJH5V3eV/3gN8CNhwPsZEOK/XC8RY4q8Bqj3xi0h9EWkYmAYuBlKrOw5jTPkEEn9cnI3E4nbhaNXTCvhQRAL7f09VPw1DHMaYcigsLASi7eJuDVDtiV9VtwLdq3u/xpjKOVLit8TvdnYEjTEhcRL/lcTGWlWP21niN8aEZN681QB8882yMEdiKssSvzEmJH/84x8BWLBgQZgjMZVlid8YE6JoIId777033IGYSrLEb4wJSe/e/YFXuemmm8IdiqkkS/zGmJCcfXY3YmKi6NixY7hDMZVkid8YExKPx0dUlC/cYZgTwBK/MSYkXi+W+GsIS/zGmJA4JX4NdxjmBLDEb4wJiZX4aw5L/MaYkHi9SrTdtFsjWOI3xoTE61Ur8dcQlviNMSHxeJToaEv8NYElfmNMSAoLhSjLGDVCOPrjN8a4kFPVE+4ozIlgh9EYE5LCQru4W1NY4jfGhMTns6qemsKqeowxISkstKqemsIOozEmJD4flvhrCDuMxpiQWOKvOewwGmNC4vOBiPXVUxNYHb8xJiQ+n1irnhrCSvzGmJA4VT1W4q8JLPEbY0LiVPWEOwpzIljiN8aExOfDqnpqCEv8xpiQWKuemsMOozEmJHbnbs1hh9EYExIr8dccYTmMIjJURDaKyBYRmRCOGIwx5aNqdfw1RbUnfhGJBqYBlwCdgetEpHN1x2GMKR+nqsea9dQE4biBqxewRVW3AojIP4ArgHUnekcPPvgCM2bkn+jNGlMr7dnTiXbt1oY7DHMChCPxnwLsKPJ3OtD76IVEZDQwGqBt27YV29EpzTjzzO0VWtcYU9yZZ6YybtwxX1XjQuFI/CWdKx5zO6CqTgemAyQlJVXodsG77hrJXXdVZE1jjKm5wnFxNx1oU+TvU4FdYYjDGGNqpXAk/u+BM0Skg4jEAdcC88IQhzHG1ErVXtWjql4RuQv4LxANvKWqdsXIGGOqSVi6ZVbVBcCCcOzbGGNqO7sPzxhjahlL/MYYU8tY4jfGmFrGEr8xxtQyohr5Q6mJyF6gorfgNgf2ncBwqpvFH14Wf3hZ/JXTTlVbHD3TFYm/MkRkpaomhTuOirL4w8viDy+Lv2pYVY8xxtQylviNMaaWqQ2Jf3q4A6gkiz+8LP7wsvirQI2v4zfGGFNcbSjxG2OMKcISvzHG1DI1OvG7cVB3EdkmImtEJEVEVvrnNRWRhSKy2f98UrjjDBCRt0Rkj4ikFplXarwi8qD/eGwUkd+EJ+pgLCXFPllEdvo//xQRubTIaxETuz+eNiKyRETWi8haEbnHP98tn39p8bviGIhIvIisEJHV/vgf88+P/M9fVWvkA6fL5zSgIxAHrAY6hzuuEOLeBjQ/at5fgAn+6QnAs+GOs0hsFwA9gNTjxQt09h+HOkAH//GJjrDYJwP/W8KyERW7P6YEoId/uiGwyR+nWz7/0uJ3xTHAGU2wgX86FvgO6OOGz78ml/iDg7qragEQGNTdja4A3vZPvw1cGcZYilHVL4H9R80uLd4rgH+oar6q/gRswTlOYVFK7KWJqNgBVDVDVX/wTx8C1uOMae2Wz7+0+EsTafGrqmb7/4z1PxQXfP41OfGXNKh7Wf9UkUKBz0Qk2T/gPEArVc0A58sCtAxbdKEpLV63HJO7RORHf1VQ4DQ9omMXkfbAuTilTtd9/kfFDy45BiISLSIpwB5goaq64vOvyYk/pEHdI1B/Ve0BXAKMEZELwh3QCeSGY/IKcBqQCGQAz/vnR2zsItIA+DcwTlWzylq0hHlhfw8lxO+aY6CqhaqaiDN2eC8R6VrG4hETf01O/K4c1F1Vd/mf9wAf4pwK7haRBAD/857wRRiS0uKN+GOiqrv9X2Yf8DpHTsUjMnYRicVJmrNU9QP/bNd8/iXF77ZjAKCqB4ClwFBc8PnX5MTvukHdRaS+iDQMTAMXA6k4cY/0LzYSmBueCENWWrzzgGtFpI6IdADOAFaEIb5SBb6wfiNwPn+IwNhFRIA3gfWq+kKRl1zx+ZcWv1uOgYi0EJEm/um6wBBgA274/MN1Rbw6HsClOC0F0oCJ4Y4nhHg74lz1Xw2sDcQMNAM+Bzb7n5uGO9YiMc/GOR334JRobisrXmCi/3hsBC6JwNj/DqwBfsT5oiZEYuz+eAbgVBX8CKT4H5e66PMvLX5XHAPgHGCVP85U4FH//Ij//K3LBmOMqWVqclWPMcaYEljiN8aYWsYSvzHG1DKW+I0xppaxxG+MMbWMJX5jjiIizYr0DPnLUT1FxhVZTkRksYg0KmNb/xCRM6oncmNCY805jSmDiEwGslX1uRJeGwYMUdV7y1h/IHCjqo6quiiNKR8r8RtTcTfgvyvTf9f1J/6+2VNF5Br/Ml8BQ0QkJmxRGnMUS/zGVFx/INk/PRTYpardVbUr8CmAOv3NbAG6hydEY45lid+YimuqTj/y4HQxMEREnhWR81X1YJHl9gCtqz88Y0pmid+YivOKSBSAqm4CeuL8ADwtIo8WWS4eyA1DfMaUyOodjam4jTgd620RkdbAflV9V0SygZuLLHcmTqd7xkQES/zGVNwnwCCcOvxuwP8TER9Ob593AohIKyBX/SMyGRMJrDmnMRXk7zf+HVW9qIxl7gWyVPXN6ovMmLJZHb8xFeQvxb9e1g1cwAGODLxtTESwEr8xxtQyVuI3xphaxhK/McbUMpb4jTGmlrHEb4wxtYwlfmOMqWX+P+tDbHXG8nQ0AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots()\n",
    "ax.set_title(\"airspeed\")\n",
    "ax.set_xlabel(\"T (s)\")\n",
    "ax.set_ylabel(\"speed (m/s)\")\n",
    "figure1 = ax.plot(ref_time[a:b], indicated_airspeed[a:b],label='indicated_airspeed', c='black')\n",
    "figure2 = ax.plot(ref_time[a:b], true_airspeed[a:b],label='true_airspeed', c='blue',lw=0.7)\n",
    "plt.legend(loc='upper left')\n",
    "\n",
    "plt.savefig('D:/taking_off_airspeed.png',dpi=1000,bbox_inches = 'tight')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "c=116052\n",
    "d=132867"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdd3gUVffA8e9JCCn0ErrSxEKvSlcEBAUVu4iI+voiqCgqKooi+EMBK/oqKBZQQWlix4IKKkgRpIUmUqSFFggkpGfP74/ZhARSNiS7m3I+z7NPZqeeHZYzd+/cuVdUFWOMMSVHgL8DMMYY41uW+I0xpoSxxG+MMSWMJX5jjClhLPEbY0wJY4nfGGNKGEv8pkQSkQEi8qO/48iNiKiInOfvOEzxItaO35jCS0QUaKSq//g7FlN8WInfmNOISCl/x2CMN1niN8WaiIwUke0iEiMim0TkOvf8O0VkSYb1VETuF5FtwDZxvCYih0TkuIisF5Gm7nWni8jbIrLQvd9fRaRuhn1d6F52VES2isjNGZYFi8jLIrJbRA669xOaYfljIhIpIvtF5G6fnCRT4ljiN8XddqALUAEYC8wQkZrZrNsPuARoDFwBdAXOByoCtwBRGdYdAPwfUBVYC8wEEJEywELgE6Aa0B+YLCJN3NtNdO+zJXAeUBsY7d62NzAC6Ak0Anrk65Mbkw1L/KZYU9W5qrpfVV2qOhvYBlyczerjVfWoqsYDyUA54EKce2GbVTUyw7rfqupvqpoIjAI6iMg5QF9gl6pOU9UUVf0L+Ay4UUQE+C/wsPs4McALwK3ufd4MTFPVCFU9CYwpyHNhTBqryzTFmojcATwC1HPPKotTSk/NYvU9aROq+ouIvAm8BZwrIp8DI1T1RBbrxorIUaAWUBe4RESiM+y3FPAxEA6EAauda4ATIhDonq4FrM6w3b95+rDGeMhK/KbYcte7vws8AFRR1YpABE6yzUqmJm6q+oaqtgGa4FTPPJZh8TkZjlMWqAzsx7kg/KqqFTO8yqrqUOAIEA80ybCsgqqWde8qMuN+gXPP7pMbkzNL/KY4K4OTzA8DiMhdQFNPNhSRdiJyiYgEASeBBDL/SrhKRDqLSGmcuv4VqroH+AY4X0QGikiQ+9VORC5SVRfOheg1EanmPk5tEenl3ucc4E4RaSwiYcCz+T0BxmTFEr8ptlR1E/AKsAw4CDQDlnq4eXmcJH0Mp8olCng5w/JPcBLzUaANzs1e3PX2V+DU2+8HDuDc0A12b/cE8A+wXEROAD8BF7i3/Q6YBPziXueXPH5kYzxiD3AZk0ciMh3Yq6pP+zsWY86GlfiNMaaEscRvjDEljFX1GGNMCWMlfmOMKWGKxANcVatW1Xr16vk7DGOMKVJWr159RFXDT59fJBJ/vXr1WLVqlb/DMMaYIkVEsnz626p6jDGmhLHEb4wxJYwlfmOMKWGKRB2/MSZrycnJ7N27l4SEBH+HYvwoJCSEOnXqEBQU5NH6lviNKcL27t1LuXLlqFevHhm6ejYliKoSFRXF3r17qV+/vkfbWFWPMUVYQkICVapUsaRfgokIVapUydOvPkv8xhRxlvRNXr8DlvgLqeRkf0dgjCmuLPEXQuvXQ+nSsGiRvyMxxhRHlvgLoeuvh0mT4OOP/R2JMbnr2LFjntZfvHgxffv2BeCrr75iwoQJZ3XcF154Ic/bTJ8+nQceeCDP2+UnzoK0a9cumjb1aBC5HFniL2RUoVEjeOgh2LcPNm70d0TG5OyPP/44622vueYaRo4ceVbbnk3iP1vZxZmSkuKzGAqSNecsZLZtg/PPd6YnTICJE+Gjj/wbkykahg8fztq1awt0ny1btmTSpEk5rlO2bFliY2NZvHgxY8aMoWrVqkRERNCmTRtmzJiBiPD9998zfPhwqlatSuvWrdO3nT59OqtWreLNN9/k4MGDDBkyhB07dgAwZcoUOnbsSL9+/dizZw8JCQk89NBDDB48mJEjRxIfH0/Lli1p0qQJM2fOZMaMGbzxxhskJSVxySWXMHnyZAIDA5k2bRrjx4+nZs2anH/++QQHB2f3Ufj6668ZN24cSUlJVKlShZkzZ1K9evVMcd55551UrlyZNWvW0Lp1a6655hoeeughwLnJ+ttvv7F69WpGjx5NlSpV2Lp1K127dmXy5MkEBATw448/8uyzz5KYmEjDhg2ZNm0aZcuWZfXq1TzyyCPExsZStWpVpk+fTs2aNVm9ejV33303YWFhdO7cuQD+Va3EX+gsXw7t2zvTrVrBoUPOrwBjioI1a9YwadIkNm3axI4dO1i6dCkJCQn897//5euvv+b333/nwIEDWW774IMPcumll7Ju3Tr++usvmjRpAsAHH3zA6tWrWbVqFW+88QZRUVFMmDCB0NBQ1q5dy8yZM9m8eTOzZ89m6dKlrF27lsDAQGbOnElkZCTPPvssS5cuZeHChWzatCnH+Dt37szy5ctZs2YNt956Ky+++GKW6/3999/89NNPvPLKK7z88su89dZbrF27lt9//53Q0FAAVq5cySuvvMKGDRvYvn078+fP58iRI4wbN46ffvqJv/76i7Zt2/Lqq6+SnJzMsGHDmDdvXnqiHzVqFAB33XUXb7zxBsuWLTvbf5YzWIm/kFm2DB5//NT788/P/CvAmOzkVjL3hYsvvpg6deoAzq+FXbt2UbZsWerXr0+jRo0AuP3225k6deoZ2/7yyy985P55GxgYSIUKFQB44403+PzzzwHYs2cP27Zto0qVKpm2/fnnn1m9ejXt2rUDID4+nmrVqrFixQouu+wywsOdnolvueUW/v7772zj37t3L7fccguRkZEkJSVl+0DUTTfdRGBgIACdOnXikUceYcCAAVx//fXpn//iiy+mQYMGAPTv358lS5YQEhLCpk2b6NSpEwBJSUl06NCBrVu3EhERQc+ePQFITU2lZs2aHD9+nOjoaC699FIABg4cyHfffZdt/J6yxF/I7NwJNWokMGzYY/Tt25eOHXuxdKklflM0ZKxGCQwMTK8DP9tnDRYvXsxPP/3EsmXLCAsL47LLLsvyQSVVZdCgQYwfPz7T/C+++CJPxx42bBiPPPII11xzTXrVVVbKlCmTPj1y5Ej69OnDggULaN++PT/99BNw5mcWEVSVnj178umnn2ZatmHDBpo0aXJGqT46Otorz2lYVU8hcvIkhIUpDRrU58033+Suu+6iUyfl7rshOtrf0Rlzdi688EJ27tzJ9u3bAc5Iemm6d+/OlClTAKfEe+LECY4fP06lSpUICwtjy5YtLF++PH39oKAgkt0PvHTv3p158+Zx6NAhAI4ePcq///7LJZdcwuLFi4mKiiI5OZm5c+fmGOvx48epXbs2AB9++KFHn2/79u00a9aMJ554grZt27JlyxbAqerZuXMnLpeL2bNn07lzZ9q3b8/SpUv5559/AIiLi+Pvv//mggsu4PDhw+mJPzk5mY0bN1KxYkUqVKjAkiVLAJg5c6ZHMeXGEn8hUq0aVK58OL0ONDIykm3bFvH22/D9934OzpizFBISwtSpU+nTpw+dO3embt26Wa73+uuvs2jRIpo1a0abNm3YuHEjvXv3JiUlhebNm/PMM8/QPu0GGDB48GCaN2/OgAEDaNy4MePGjeOKK66gefPm9OzZk8jISGrWrMmYMWPo0KEDPXr0yHRjOStjxozhpptuokuXLlStWtWjzzdp0iSaNm1KixYtCA0N5corrwSgQ4cOjBw5kqZNm1K/fn2uu+46wsPDmT59Ov3796d58+a0b9+eLVu2ULp0aebNm8cTTzxBixYtaNmyZXprqWnTpnH//ffToUOH9PsH+aaqhf7Vpk0bLe6OHFG99VbV1157SwHdvHmz1qxZU9u1a6cHDqhef72/IzSF0aZNm/wdgsnCokWLtE+fPj49ZlbfBWCVZpFTrcRfSKxeDRdfDH/8sZhzzjmHCy64gBEjRvDnn39y+HAE4eGwYYO/ozTGFAeW+P0kNtZppqkKv//uJP7WrZXff/+dLl26ICLcdtttAHz77bfccAN88YWfgzammHj++edp2bJlptfzzz9fYPu/7LLL+OabbwpsfwXNWvX4SaVKsHgxBAVB165w9dVw7bX7OHDgQPoj8DVq1KBRo0YsW7aMESPgpZf8G7MxxcWoUaPS28mXRFbi9zFVOHAA2rVzumW47z544AHn5u3WrX8C0LZt2/T1W7ZsyZdffsmSJb9SpgzYQEvGmPyyxO9DqvDuu1CzJjz4IHz4oVOKf/JJmDbtBNdffz0AzZs3T99m6NChAIwbN47GjSGXBw+NMSZXVtXjA2PHOm30AwNh6lQYOBB694aKFcH9VDrHjjlPLHbt2jVTk61u3brx9NNP8/zzz3PddYdZsyacXFqkGWNMjrxW4heREBFZKSLrRGSjiIx1z68sIgtFZJv7byVvxeAPLpfTOift6fmUFPjpJ6hTx+l0bd8+p9O1ihUzb7d8+XLq1KnD4sWLz9jnPffcg4iwceN0Vq70/mcwxhRv3qzqSQQuV9UWQEugt4i0B0YCP6tqI+Bn9/tio18/aN3aqbNfsgQWLoS+fZ2qHZcLQkIyr+9yuVi9ejURERE0a9Ysy8ez69atyxVXXME337zFrl3WY5spPKKjo5k8ebJfjl3c+sj3Ja8lfvfzA7Hut0HulwLXAmnPQn8I9PNWDL525AgcPQpvvQVxcdClC/TvDwMGOMvTcnpMTAzNmjVj2rRpvP3227Rt25Z169al90aYlTvuuIPdu/8lJeUwUVE++DDGeCC7xJ+amur1Yxe3PvJ9yas3d0UkUETWAoeAhaq6AqiuqpEA7r/Vstl2sIisEpFVhw8f9maYBcLlghEjYPBgpy5/zhyIiYEdO5xqnowWLVpEREQEd999N/fff3/6/JwSf+/evQEIC1tLAfbOaky+jBw5ku3bt9OyZUvatWtHt27duO2222jWrNkZJeGXX345vdOz7du307t3b9q0aUOXLl3S+7fJytdff80ll1xCq1at6NGjBwcPHgQyj6Z155138sgjj9CtWzeeeOIJfv311/T2+a1atSImJobFixfTtWtXrrvuOho3bsyQIUNwuVwA/Pjjj3To0IHWrVtz0003ERvrlFlXr17NpZdeSps2bejVqxeRkZHp81u0aEGHDh146623Cvy8eptXb+6qairQUkQqAp+LiMe/h1R1KjAVoG3btoW6fmPxYhg/Hnr2dG7cAtSokfW6H374IXfeeWemebVr12bfvn306NEj22NUqlSJCy64gBMnFrBkyRW4R64zJpPBg2H//oLbX61aToOE7EyYMIGIiAjWrl3L4sWL6dOnDxEREdSvX59du3blEOdg3n77bRo1asSKFSu47777+OWXX7JcN62PfBHhvffe48UXX+SVV145Y720PvIDAwO5+uqreeutt+jUqROxsbGEuOtYV65cyaZNm6hbty69e/dm/vz5XHbZZel95JcpU4aJEyfy6quv8uSTTzJs2DC+/PJLwsPDmT17NqNGjeKDDz7grrvu4n//+x+XXnopjz32WJ7OaWHgk1Y9qhotIouB3sBBEampqpEiUhPn10CRoOrcrA0Kyjz/hhuchD9iRNbbffzxxyQnJ7Nhw4b0PtMvv/xyJk+ezJEjR7jwwguJjIxM78c7O+3bt2fBgtmEhb0GFHxXraboyylJ+8LFF1+cbR/2aWJjY/njjz+46aab0uclJiZmu35J6SPfl7yW+EUkHEh2J/1QoAcwEfgKGARMcP/90lsxFLQvvnCaZmYc3W7+fOjT51QrHlXlnXfeQUS499572bZtG3fccUem/QwaNIhnnnmGhg0bcsEFFwCcMbBEVi655BI+/PBDVOM4dqwMlYpVeyhTHGTsp75UqVLpVSlAej/6LpeLihUrejxMZEnpI9+XvFnHXxNYJCLrgT9x6vi/wUn4PUVkG9DT/b7Q+/df+N//nFY5R4868+bOdZpoPv30qfWefPJJhg4dypAhQ9izZ0/6lymtS9YdO3Ywffp0GjZsmOcY0rqkbdhwPf/3f/n+SMbkW7ly5YiJiclyWfXq1Tl06BBRUVEkJiam911Tvnx56tevn943vqqybt26bI9RUvrI9yWvlfhVdT3QKov5UUB3bx3XW+69F+rWhZtvdkr3Dz0E99wDLVueGh1ry5YtTJw4MX2bc889F4CePXsyf/78fMfQrFkzQkNDKVVqDps3d8j3/ozJrypVqtCpUyeaNm1KaGgo1atXT18WFBTE6NGjueSSS6hfvz4XXnhh+rKZM2cydOhQxo0bR3JyMrfeeistWrTI8hhpfeTXrl2b9u3bs3PnzlzjmjRpEosWLSIwMJDGjRtz5ZVXsmzZsvQ+8jds2JB+ozcgICC9j/y0Kqdx48Zx/vnnM2/ePB588EGOHz9OSkoKw4cPp0mTJkybNi19APRevXrl8yz6QVZ9NRe2l7/749+9W/W225xpl0u1fn3Vc85RXbLk1Doul0v79eungP79999aoUIFxWm+qgsWLCiwWLp06aLt27fXq69WTUoqsN2aIsr64/ecP/rI9yXrj7+ALVnidLEATlv8775z6vvd94KIi4vj5Zdf5osvvuDmm2+mUaNGbN68mUOHDpGUlJQ+Ik9BaN++PX/99RcXXJDCxo0FtltjTAliffV4YP16uPXWU+/d92MBiIqKol69esTGxnL99dfz0UdOnzs1a9b0Sizdu3fnpZdeIjBwLatXt6VlS68cxhife/75588YE/emm24qsO6TL7vsMi677LIC2VdRZ4nfAxs2wJgxsHjxYhYuXMi4cePS7+p/+umnxMbGcu211zJz5kyCg4O9Gku3bt2oWLEiERFfU7p029w3MMWeqhb5ViZgfeTnh1Or4zmr6snFL7847feDg+Guu+7ihRdeYMWKFQBs2rSJl156iVq1avH555+nPyTiTaVLl+b666/nl18+ZteuJK8fzxRuISEhREVF5fk/vik+VJWoqKg85R8r8ediyBD48EMXAwYMTH8S8Z133iE4OJjW7v6Rx48f79MS14ABA/jgg2ls2XIMqJ7r+qb4qlOnDnv37qUodGtivCckJCTXB0AzssSfg9hYaNwYTpxYyCeffELdunXp0KED06dPZ9asWYDTv87DDz/s07icC44SExOLJf6SLSgoKNcnZY05nVX15GDRIrj8cnjvvfcIDw9n69atTJgwgdKlS5OQkMBnn31GRESE1+v1T1exYkVq1KiR3pGUMcbkhSX+HMyZA1dfrSxdupQrrriC4OBg6taty/r161m3bl36UIn+cNFFF3Hy5FEbg9cYk2eW+LPx77+QnAywi8jISDp27Ji+7IILLsg0Lq4/OIk/gm3b7KaeMSZvLPFnY/JkZ9Ssn3/+GYAuXbr4OaLM6tevT1LSUpYssSK/MSZvLPEDUVGwd++p99dcA19/De3apTB69GjCw8NzHCTFH2rVqgWsY8WKeH+HYowpYizxA6+8At3d3cZ98w2ULQsbN8LChT8SGRnJCy+8QEBA4TpVTuLfxtatrlzXNcaYjApXNvOT7dudh7R27IB334V33nH65Jk7dy6VK1dmYNqwWoWI001tKtu3l+bECX9HY4wpSkp84leF48fhq6+gYUNnhK1y5ZxBI7755ht69Ojh8+aanqhfvz7BwcE0aPAnO3b4OxpjTFFS4hP/7t1OP/sXXghHjjhVPQCffPIJR44cYciQIf4NMBulSpWicePGnDixkT17/B2NMaYoKfGJf+nSU90rV6niVPGoKq+//jrNmzcv1L35NW/enMjIleze7e9IjDFFSYlN/DExkJoKr70GzZsf4t9//wWc4dX69OnD+vXrGT58eKHu9bBFixZER69jy5aT/g7F+MH48bBiBbz0EixfDvHx8MMPTvWlMTkpkYl/61YoXx6aN4f//AcGDepJ27ZtSU1N5bnnnuO7777jmmuuOWOQ9MKmZcuWwG7WrDni71CMD6nCq6/CU09B375O4aVDB5g3zxkwaPx4f0doCjspCt25tm3bVletWlVg+5s0CRo1gpo14fzzYylXrhwAEydO5JlnnuHaa69l9uzZhbq0D86vk/Lly1Ojxip27ixczxkY70m7L/X6607B5ZdfnFJ/dDRcdRWsWweffgr798O55zrNk03JJCKrVfWMgTtKZIn/11/h0kuhdWtYuXJl+vwnnniCpKQkHn/88UKf9OFUz4zx8fb0ri988IHT9PfgQf8cPybGSeibN0OTJnDjjVCmDFx9Nfz2G3TuDAMHQr9+8PPPzjrz5vknVlO4eS3xi8g5IrJIRDaLyEYRecg9f4yI7BORte7XVd6KISvbtzvdLaeVgv744w8Afv31V9q1a8eNN95ImzZtfBlSvpxzzjkkJR3npFXze9WBA/Dyy04yrVEDEhN9H8PHH8NttzlVlS+/DLVqZV4+ebITX8+e8MQT0LYtTJ1qdf7mTN4s8acAj6rqRUB74H4Raexe9pqqtnS/FngxhkzWrHH+47z00ql5y5Yto3HjxnTt2pWVK1cyd+7cIlHaT1OnTh2Sk1ezZo2/Iyl+TpyATZucxPnTTzByJCxcCIMGOfPB6cjv2DHvx7Jvn5P4e/RwEn9CwjoaNmxIWFgYIkKHDh346KOP6NKlCwcPLmP2bJgxw2mm/Ndf3o/PFC1eG4hFVSOBSPd0jIhsBmp763i5SUqCRx6BBQugcmUlImIjZcqU4bfffqN///7+CivfGjRoQGzsO7z55sN07mzj6hSkQYOcp7mrVXNaz2zd6twX2rkTnnkGLrvMSfq//ALLlnk3lnvucUryq1fDd9/FM3lya+BUdx3Lly9n+fLlAHTs2JETJ05Qrlw5brnFqdosQj9ijQ/4pI5fROoBrYAV7lkPiMh6EflARCpls81gEVklIqsKYli5P/5wSktVqsCECRNo1qwZDRo0ICgoyOcjaBWkBg0aANtJTIxl61ZnXkICuKwLn3w5eNB5gnv/fqe0HxPjJH1wnvD+9luYPh3+/BNCQpxqlzffdJpUFjRVp+nxc89B6dLKwYMnSEv6H330EVdddWZt6bvvvsstt9zCjh3zWLu24GMyRZvXE7+IlAU+A4ar6glgCtAQaInzi+CVrLZT1amq2lZV24aHh+c7jh9+gF69nOmPPvoIcKpJPvvsMy666KJ8799fnMQPHTtuZMYMZ17fvtC+vXMvI+1nfmioM8ZARrfeCqNHw6FDTmJLk5BAev8/JfUCsno1tGsHjz/unLdff91OmTJlaNy4MR98cC9btpygShUn6SckOHXrw4bBZ5/ByZMQF1dwsWzbBuefD1OmTGHBgkWULbsLgNdff52BAwfyySefsHTpUqpWrcrYsWNp164djz76KHPmzOGxx+7kiLX2Nafxat2AiAThJP2ZqjofQFUPZlj+LvCNN2NIs2YNPP88HD58mC1btjBx4kQef/xxXxzaq9ISf6lSq/nzz05ERUHVqk6zvksucbqb3r4dzjvPubcxaBC0auUk+9RU55fQsmVO6fbGG51qjXvvherVnRuEDz/sVHd07+60EKmU4fdZVJTzC6o4ioiAiy92qnMApk6dRlxcHJs3b2bz5s20bt2cPn3up2JF575RaKhzsbz5Zudp8H/+gVmzYO5c59dA377gSQevJ086LXXSrFvn/OK4/PJUbrjhPmAahw6d5JprruHBBx8EoEKFCnTs2JFDhw4hIoSGhvLnn3+693eS6OhDJCRUIySkYM+RKbq82apHgPeBzar6aob5NTOsdh0Q4a0Y0uzbBxUrgmoqr732GkCmEbWKsqpVq1KzZk2WLv2NypWdpH/ppc79jI0bnVL9d9/B7bc7NyQvvhg++ghq13bagIeEQFCQ01Jl8mTo398p7U6Y4DQLrFzZuXCsXAmPPuo0JXz3XeeCUK+e08FdcrIzfkFR8vrrTjXN6Y4fd6p0Zs+Gr7+eSO/evZk1axYvvvgiHTp0SK8W/OOPP7jttr0EBr7Pc889TmCgc1G84gp4+20ngd9wg/N37FgYPhymTHHOU/Pm8MknMGKEcz4Bxo2DK690fmXMn+/M27PHuZCMGAGvvHKNO8KhwMNZjg+R1ijhvvvuY+zYsTz77LNUqlSJo0d/45NPnJu9qakFfCJN0aSqXnkBnQEF1gNr3a+rgI+BDe75XwE1c9tXmzZtND/69VPduFF1+PDhCmjv3r01OTk5X/ssTO677z4NCgrSv/7ar7/9pnry5Kllb7yh2rSp6rffqq5bp/rSS6qgun27s/zWW1U7dXKmk5JUjx49tW18vOo776h266Y6fbqzn7Ztne0HDHD22aeP6oUXOvM2b/Ys3tRU1fffV3W5Ms+PilJNScl6m6ZNVZct82z/uTl6VLVdO9Urr3Q+/+uvO7Hv3Kl6883O98XdCDLTa/To0aqq2qpVqzOWjRkzRnfv3q1btrj06aedz3Lw4Kljtm+vOnGi6pgxqn/9pVqnjmqPHs65XbXKiUVVNSHBmX7oISeGZctUf/jh2BnHe//99z36rAMHDtSKFe/RUqVUr79e9fHHVU+ccP5tL7pINTq6YM6pKZyAVZpVfs5qZmF75Sfxnzypeu21qn/99ZcCOnjwYHWdnnGKuE2bNmlAQIAOGDBAVVWff/55nTp1qqampmpsrPOv/PDDr2lUVJTGxKg+/LCqy+XSL774Qnv2PK7dumW/7wMHnO0PHlQ9dky1Z0/V3393kqTLpTpvnnMR+fNP1XvucbaZN0/1wQdP7ePff1UTE1VXrlRt3Vo1PNzZ57ZtmY/Vp4/q1Vc7F4WMVq1S7d9ftWvX7C8Mqs4FpWtX1a1bnfeJiaodO6qOHKk6aZLzXXC5nP1/9JHq00+rtmihWrq0c2wn2TvbP/bYm2ck219//VVVVe+4444zlqW9xo4dm+u/l8t16kIZEaFavbrqqFGnlnfs6CxfssR5//nnnyugP/74Y/pxtp1+8rLx+uuvK5TSiIj9um+fs9+HHlJ9803Vm25Sfe45j3ZjiqgSmfj37lX9739VJ0xQfeqppzQwMFAPHz58Vvsq7O666y4tW7as/vrrr+nJQUS0efPm2rTpRgW0U6dOmuLOnKeSyDQtX/47Ha2r5oYAACAASURBVDFihCYmJuquXbv06quv1rvvvltXrlypY8eO1d9+c+mMGTM0Kioqxxh69XJ+CXTv7iTwd95x3pct6yTeG25QjYx0knP//qqzZqm2auWUjv/5R/Wuu1R37FBt0sQpiQ4d6ux35Egn+T/zjOqCBapPPOGUnFNTTx27e3fn23zLLc7yw4edYw4f7sxv08aJ45FHnPfR0c72J0+qxsQ4+/jhh1PTN998s9atW1dvuOEGfeWVV/TTTz9NLzDs27dP586dm2XiL1OmjI4fP17/85//qMvl0n379uVa0ADVH3889b56ddWnnnKmXS6XPvrooxocHKwJCQnpx/HU0qVLFdAvv/xSVVVjY1XPPVe1Xj3nF97ll5/5y8sUHyUy8aemqlaurLphg+qFF16ol19++Vntpyj48ssvMyWg5s2bp09Xq1Yt07Jx48ZleF9OoWy2JVhA33nnHQX0zjvv1AMHDuj69ev1xRdf1J49e+qhQ4dUVTU+Pl5dLqfqZNs21T/+cH4dTJ6s+vffqs2anapSUnWqnR57zPkG3n676vnnqy5d6iwbN061cWPVatVU169Xveoq599yyRLnQlGz5qkqp8aNVTdtci4q0dFOKb9GDdXQUNVGjZwLwNq1qo8+eqpEf/Cg6vLly3Xnzp1nnMfdu3fr0KFDtXr16nrdddfleM6vvvpq7d+/v27YsEH379+vX331Vabz1rBhw/Tpxx57TP/++289duyYTpw4UX/88UddsGCBRkZG6v/+p/rNNz/pE088od26ddOuXQ/oZ5+p7t27Vxs0aKCAXnrppaqq2qxZMw0PD/f4e3Hy5EkNCAjQoUOH6nZ3/d727U7SV1UdPNg5l3FxHu/SFCElMvGrOqWZTZs2KaBvvvnmWe+nsIuJiUlPMk2aNFGXy6UpKSkaHBysgNauXfuMhH7DDTekT5crVy59eu7cufr000/neDFIe/Xr109vvvlmBfT333/PNr7PP3fuN6iqpqam6vr1m7R2bZfefbdzcdi69VTJ0+VyEv3y5U7Sv/LKVB09erQuWxah4CTypCTnV8CoUaqBgU6VSUanVwm99dapxL969WoFtGrVqjp27Fj94osvVFU1Ojo6/bOk1dvnxOVynVGar169eqbzU6NGjRzP37nnnqs//PBDpnm1atXSTz/9NNO8119/XVVV//33X42MjMzl25BZxpiefPJJdblcGhUVpWPGjNH9+w/rpZc65yWnajRTNJXYxK/qVIMAunfv3nztp7BbtGiRPvXUU5lKsh06dFBAb731Vv3ss89URBTQV199NVPVQXR0tB49ejRT3fHevXu1Xbt2Z/yCSHtlTJKAdu7cWVVVly1bph999FH6ftKql3bv3p3+bwHoFVdE6A8/nFqnZ8+eOnz4cF2W4S4uqN5xx1oFNCwsTGvXTs2UoGJinGqitH1MmDBBN27cqKqqSUlJ+r///U/Hjx+vzz//uw4a5NzYfeCBBzy6qH311Vd5/jeYPXu23n333bp161Zdv369qjpVQ9WqVdPevXvrAw88oJ06dcryeLt379bZs2dnmjds2DB97733ND4+Ps+xpBk6dGimfQ4cOFD/85//KKA33nij9u/vnOfTL56m6CuxiX/9+vUK6P3333/W+yjKpk6dqoAOGTJEVZ2b3LfccoseP35cVTU9GWRn8uTJCugbb7yhv/32m3bu3FnHjh2rixYt0qioKJ04caLOmjVLH3zwQQX0+uuvT9/nyy+/rMOGDdNKlSrp008/rXXr1s2UgG6//XZ99NFHddu2bWdUVa1bt06nTZum06alaM+eV6fPHzlypM6ZM0e/+eab9AuKy+XSQYMGpa/TpUsXdblcOnbs2Ez7PP/887Vfv34KaN++fc+IJ+MxFi1aVKCNAFJSUtL398svvyig48ePT29pVrFiRVVV3bJlS3ocI0aM0ISEhHwfOyYmRvfu3atz5szJ8vNu2bJLJ01y7p+Y4qVEJv6EhAS9/vrrVUT0wIEDZ7WPou7o0aPar18//fvvv7Nc3q5dOz3vvPOy3T4pKUlnzJiRawI6duxY+q+JnF5vvfWWNm7cOMtlVapUOaOKY8aMGRoWFqYPPPCANm3a9IxtOnfurF9//bVHJfiMr4ULF+qhQ4f0+PHjeujQIT1x4oS+9957um3bNp+0+kq7aK1YsUIB7dWrl6rqWd3AzYuIiIj0/S9evFgBLV++vH7+eYq+/bZXDmn8qEQm/oEDByo4N9lM1iIjI3NtreOpw4cPa9++ffXrr7/WGTNmaK9evXT8+PE6f/789GSzb98+VVV97rnnzkjGEyZMUFXVq6666oxl8+fP1//+97/ZJvKyZcvqnj17dOvWrZna2S9ZskQXLVqUad1PPvmkQD5vQUhJSdHnnntOd+/enT6vT58+Om7cOK8dc926dRoXF6culyvDeWmiXbrEaDFt9FZilcjEn3YTb1TGRtLGL9566y2dOHFi+vuMNzRff/11XbNmzRkl7REjRqSvExMTo1OmTFFABw0apLGxsbpy5Urt3r27gtPGPaPT95V2kxvQEydOeO+DFjEHDx50n5cgBad9vyk+SmTiV1Xdv39/+s9qU3ikpqbqqFGjMt0EPl1atURaM9wDBw5o9+7dM1VbHTt2LL2ZYk5mzZql9957ryX9LFx55ZUKaI8ez+qTT/o7GlOQskv8JXLMXVM0qCpz5syhZ8+eVK5c2d/hFFuqyuWXX058fCnq1l3I7Nn+jsgUFBtz1xQ5IsItt9xiSd/LRIS6deuyYsVPHDxYgP1Jm0LLEr8xhkaNGgHwxx+/+zkS4wuW+I0xXH311QAkJ8d4ZRQxU7hY4jfG0Lx5c5555hlgPzt3Jvk7HONllviNMQA0bNgQ2MuaNYf8HYrxMkv8xhgAzjvvPGAf69bZIL3FnSV+YwwA5557LrCP7dsT/R2K8bJcB1sXkWpAJ6AWEI8zRu4qVXV5OTZjjA/VqFED2MfevYX/2R6TP9kmfhHpBowEKgNrgENACNAPaCgi84BXVPWELwI1xnhXUFAQ1aolc/hwkL9DMV6WU4n/KuC/qrr79AUiUgroC/QEPvNSbMYYH6tduzJ799qP+eIu28Svqo/lsCwF+MIrERlj/KZWrVps357g7zCMl+V6c1dEHhKR8uJ4X0T+EpErPNjuHBFZJCKbRWSjiDzknl9ZRBaKyDb330oF8UGMMflXu3ZtEhLiSE31dyTGmzxp1XO3ux7/CiAcuAuY4MF2KcCjqnoR0B64X0Qa49w3+FlVGwE/u98bYwqBWrVqkZS00x7iKuY8Sfzi/nsVME1V12WYly1VjVTVv9zTMcBmoDZwLfChe7UPcW4WG2MKgVq1agHbWbXqmL9DMV7kSeJfLSI/4iT+H0SkHJCnuz8iUg9oBawAqqtqJDgXB6BaNtsMFpFVIrLq8OHDeTmcMeYs1a5dG9jDpk3H/R2K8aJsE7+75Q7Af3CqY9qpahxQGqe6xyMiUhan5c/wvDT9VNWpqtpWVduGh4d7upkxJh/q1asHHGD79lh/h2K8KKfmnMtFZC/wPfC9qkYDqGoUEOXJzkUkCCfpz1TV+e7ZB0WkpqpGikhNnOcDjDGFQFri373b6viLs2xL/O5RWx5yv50kIn+KyGsicoWIBOe2YxER4H1gs6q+mmHRV8Ag9/Qg4MuzC90YU9DCwsKoWjWVyMhcb+OZIizHOn5V/VdV31bVfkBH4GugB/C7iHyby747AQOBy0Vkrft1FU6LoJ4isg3nATBPWggZY3ykfv3KREf7OwrjTbn21ZNGVZOBX9wvRKR2LusvIfvWP909Pa4xxrfq1avLxo02Gktx5skDXH1FZI2IHBOREyISIyInVHWfLwI0xvhW7dq1SUyMIyXFOmsrrjxpzjkJpy6+sqqWV9Vyqlrey3EZY/ykdu3apKbuYssWa9lTXHmS+PcAEapql39jSgDnIa5INm70qPGeKYI8qeN/HFggIr8C6SM0nNZSxxhTTJxzzjnAr2zadBSo5+dojDd4kvifB2Jx+uIv7d1wjDH+1qhRI+AgW7da057iypPEX1lVc+2N0xhTPFSvXp2wsFh27ozzdyjGSzyp4//Jk26YjTHFg4jQsGFZIiOtb+biypPEfz/wvYjEZ2zO6e3AjDH+U79+GNHRIf4Ow3hJronf3XwzQFVDrTmnMSVD7doVSUgIYO1af0divCGn3jnr5bShe0SuOgUdkDHG/6pXr05yciitWvk7EuMNOd3cfUlEAnA6UVsNHMZp2XMe0A2n24Vngb3eDtIY41vVqlUDKvs7DOMlOQ22fpN7qMQBwN1ATSAOZyStBcDzqmqjMhtTDFWvXh2o4O8wjJfk2JxTVTcBo3wUizGmkHAS/0zgcXbs2EGDBg38HZIpQB73zmmMKTmcxP8E0Jhjx2r6OxxTwDxpzmmMKWGcxA+QQmys1egWN5b4jTFnKFu2rHsqhiNHEnNc1xQ92Vb1iEjrnDZU1b8KPhxjTGEgIowZM4YxY2KIirLxd4ubnOr4X3H/DQHaAutwRtRqDqwAOns3NGOMPw0YMIAxYz6zEn8xlNNg691UtRvwL9BaVduqahugFfCPrwI0xvhHhQoVgOMcPpzs71BMAfOkjv9CVd2Q9kZVI4CW3gvJGFMYOIn/IAcO2BhMxY0niX+ziLwnIpeJyKUi8i7OQ1w5EpEPROSQiERkmDdGRPaJyFr366r8BG+M8Z7SpUsTFHSIQ4eC/B2KKWCeJP67gI3AQ8BwYJN7Xm6mA72zmP+aqrZ0vxZ4GqgxxvcqVjzOwYOh/g7DFLBcH+BS1QQReRtYoKpbPd2xqv6WW0dvxpjCrUaNVKKirHvm4ibXEr+IXAOsBb53v28pIl/l45gPiMh6d1VQpXzsxxjjZTVrluPkSXvcp7jx5F/0WeBiIBpAVddy9iMwTwEa4twcjuRUk9EziMhgEVklIqsOHz58loczxuRHjRrVSUy0dvzFjSeJP0VVjxfEwVT1oKqmqqoLeBfngpLdulPdTUjbhoeHF8ThjTF5VL16dRITE1G1lj3FiSeJP0JEbgMCRaSRiPwP+ONsDiYiGXt7ug6IyG5dY4z/Va9eHVUXx48XSNnPFBKeJP5hQBMgEfgEOI7TuidHIvIpsAy4QET2ish/gBdFZIOIrMcZzOXhs47cGON1TmdtyoEDB/0diilAnrTqiQNGicgLqnrS0x2rav8sZr+fl+CMMf7lJP7D7Nx5hAsvvMDf4ZgC4kmrno4isgn3Q1si0kJEJns9MmOM3zmJP5qdO6P9HYopQJ5U9bwG9AKiAFR1HdDVm0EZYwqHqlWrAsf45JMGLLDHLYsNjxroquqe02aleiEWY0whU6VKFeAYS5deRJ8+/o7GFBRPhl7cIyIdARWR0sCDeNBXjzGm6AsODqZ06ThcrqOEhFQAAv0dkikAnpT4hwD3A7WBfTgPX93vzaCMMYVH+fKppKRUJjZ2lbXnLyY8adVzBBjgg1iMMYVQxYrKkSMAx/nnn39o1KiRv0My+eRJq54GIvK1iBx2d7P8pYg08EVwxhj/q1JF3FPxHDwY5ddYTMHwpKrnE2AOUBOoBcwFPvVmUMaYwqNcuTj31DF27z7h11hMwfAk8YuqfqyqKe7XDMAq+owpIcLC0hL/UfbujfdrLKZgeNKqZ5GIjARm4ST8W4BvRaQygKoe9WJ8xhg/CwoS4AXgJPv31/Z3OKYAeJL4b3H/vfe0+XfjXAisvt+YYqxUqVLAKOBeDhywLpqLA09a9dT3RSDGmMKpXr167qmjHD7s8mcopoB40qrnJhEp555+WkTmi0gr74dmjCkMxowZw/Tp06lVK4SjRyX3DUyh58nN3WdUNUZEOuP02fMh8LZ3wzLGFBYhISEMGjSIKlWE48ftyd3iwJPEn9YvTx9giqp+CZT2XkjGmMIoPDyAmBj7r18ceJL494nIO8DNwAIRCfZwO2NMMVK1agjx8fZfvzjw5F/xZuAHoLeqRgOVgce8GpUxptCpVq0MiYmeNAQ0hZ2nI3DNz/A+Eoj0ZlDGmMInPLwCKSmBJCcnExQU5O9wTD7Y7zZjjEfKli0DQHy8Pb1b1FniN8Z4JCwsDFDi4uJyXdcUbpb4jTEecRK/WOIvBryW+EXkA3c3zhEZ5lUWkYUiss39t5K3jm+MKVhlyjhVPSdPnvRzJCa/vFninw70Pm3eSOBnVW0E/Ox+b4wpAqyqp/jwWuJX1d+A03vuvBbnyV/cf/t56/jGmILlJH4s8RcDvq7jr+5uDprWLLRadiuKyGARWSUiqw4fPuyzAI0xWXMSv4uYGEv8RV2hvbmrqlNVta2qtg0PD/d3OMaUeE7ijyM62rpmLup8nfgPikhNAPffQz4+vjHmLDmJ/yTHjlniL+p8nfi/Aga5pwcBX/r4+MaYs5RW4j9xItnfoZh88mZzzk+BZcAFIrJXRP4DTAB6isg2oKf7vTGmCEgr8R8/nprruqZw81qPS6raP5tF3b11TGOM94SGhgInOXHCf4k/NhaioqBuXb+FUCwU2pu7xpjCpVSpUgQGJhIT47/hF6dOVVq00PT3t9wCLhsNMs8s8RtjPFa6dAqxsf7LtHPnfkts7In093PmwNat0QCowoYNcOyYv6IrOizxG2M8FhycQmys5r6ilyxfvoHU1Gji4uIYNsx5vufpp18DICAAWrRQ3nsvxW/xFRWW+I0xHgsJceHfXpkrAicoUyaMN990nu+ZP38sSUlOE1NVYfTof/0YX9Fgid8Y47GQEPyW+FUVqAA0OGPZ0qUr3VNHSUho6MuwiiRL/MYYj4WFCQkJ/jl2bGwsUA4oc8ayUaMmuaeCgMXWg2guLPEbYzzmJH7xy7GdPrtOP/Y64H0OHGjmfl8OuIyBAw/6NLaixhK/McZjYWGBJCb6LvEfPgzjxzvTR48exUn8Y4F/cYb+HgNUZefOZzNstY3vv9/ssxiLIkv8xhiPlSkTQFKS79LGokX/8NRTTlPNP/9M6yNoDBAFfAN8Adx/2lY3Eh+/ljfe8FmYRY4lfmOMx8qUCSQ5OdBnxxs9+iUAZsyI4r77OgIpjBo1CvgA+J6FCxcC+4A/gF/cW8UArXnoIdJb+5jMLPEbYzxWrlwpkpO91tPLGXbujAJg0KDZ7jmHefDBB4G3gPn06NGD33//HegIXA6UplOnlsBFAIwatdBnsRYllviNMR4rVy7Ip4k/IKACAKppg78cpFKlSvzwww8sWbIEgEsuuQSYCXxNgwbn0Lp1I6AekMwHH9hYHlmxxG+M8Vi5ckG4XEG4fNRBjkgF99Rd7r/HCQoK4oorrqBTp04ABAUFAbcD17B9+3bOOSct2X/L0aMXk5iY6JNYixLfXbqNMUVemTJBQGkSExPdvXV6V6lSlYE1QCv3nKyfHpsyZYr7AS8oX76ce+4woCrz5v3LgAEDvBxp0WKJ3xjjsdDQEAASEhJ8kviddvn/BW4A6gB/ZbnWkCFD0qfLlEl7wGs/cIzg4OAst/n4Y6eb56FDCzDcIsISvzHGY06yFxJ89PhuYmIIcBR4CoALLrjAg22cqp2OHdvzxx+pnDyZdVXPyJG/ExdXlqpVW3HTTQUVcdFgdfzGGI+FhJwq8ftCcnIoQ4feztq1a1mzZg1Lly7NdRunxB9EvXr1gJMcPZp14t+/fz/R0ce4+Wb4++8dBRp3YWeJ3xjjMSfxK/E+6KktOTkZ1TKEh4fQokULWrZsSZUqVXLdrk+fPgwZcg//93//R06J3+E8hTxs2NiCCbqIsKoeY4zHfFXiP3ZM6dXLBQRToUJInrYtV64cU6ZMITU1FYglOjoo121+/PFD4uLiCQvzxX0L/7MSvzHGY75K/C+88B5//unclK1QoUIua2ctMDCQUqWSOH7cs4FZ/u//JpzVcYoiS/zGGI/5KvFPnvype6ovFStWPOv9BAcnc+JETs8cnOpwzpcPpvmbXz6piOzC6VAjFUhR1bb+iMMYkzdOq54Yryf+hIRTfe6fbYkfICQkhZiY1DPmp7X5h1PVSKGhlc76OEWNPy9x3VT1iB+Pb4zJI6fE7yIuzrudn7lcacl+DYGBZ98pXGhoKjExZ853Om8LAEoDq4C2bN1a7ayPU9RYVY8xxmNO4k8hLs7b3SBUdv99h1atWuW4Zk7CwpSsBuNyRuhyWig5TwbD3Ll1z/o4RY2/Er8CP4rIahEZnNUKIjJYRFaJyCpn5B1jjL+lJf7sHooqOGnNNmPyVcdfpgzEx2dOc+vXH6ROnTD3u2AgEegMfHnWxylq/JX4O6lqa+BK4H4R6Xr6Cqo6VVXbqmrb8HDrYc+YwsB3iT+txB+br72ULSvEx2euKnrggWeJjw/BST/BQALnn18jQ4dwxZ9fEr+q7nf/PQR8DlzsjziMMXnj9HuTQmLimTdMC4pz49VJ/AEBcTmvnIty5YSEhMyJf8uWtPF4S5GW+Dt0aEpQUMm5uevzxC8iZUSkXNo0cAUQ4es4jDF5V6pUKSCFpCTvdcs8evRxoAMAQUH5az1Uvnwgx4/ffNrcyhmmQ4FEKlUKIilpcIbWPsWbP0r81YElIrIOWAl8q6rf+yEOY0we+SLxjxtXEWgAQFBQ/qqUtmxZhWoNli1blj4vJiatt85fcW7wJrBtm9Pr56xZs/J1vKLC54lfVXeoagv3q4mqPu/rGIwxZyct8Scne6dk7JS4T5XyS5fOX+Jfu9YZpeu3335LnxcaWge4zP0KApLZsGE1AN9882O+jldUWHNOY4zHAgIC8GaJ32lff+r+QenS+buX0K1ba2BhppZBLlcF+vTpzMGDBwEXEMADD9wPvE/Zsufk63hFhSV+Y4zHRAQRl9dK/E6vn2kXlTp06HB+vvY3ZcoUIIHFiwO4885/AYiNbUC1aoHu7ptTgEAeffRRoCLR0bXydbyiwhK/MSZPAgK8l/jj4uJIq+qZOPFBPv7443ztz+nuIYFZs/bz4YcViI+PJzW1OsnJ+9zdT6QCge5fMsqcOflrPlpUWOI3xuRJQIDLa1U9Tonf6WOhR48eGYZRPDthYWFAHE5nbCF89tl2oDUNG9ZxJ3sn8TveA3bl63hFhSV+Y0yeBAbi5aoepwuvoKDc+9HPjVOqj8NJdSEMHDgfgIYNG56xbrdu7QgPr5/vYxYFlviNMXni/aqeKCDtRm/+BAUFERCQiNN6B9JSnvNLAOBS4B3AafOfkJD/i01RUHI6oDbGFAgn8Xtn306J36l+qV17X4HsMygomcTE0u53Tv/7aVVIGzfOTW/xU6FCKRITLfEbY8wZAgOVFM8Gtcozp59/ZenS36lRo0aB7LNUKSEx/XGAzCX+xo0bp69XrlwQKSln3wV0UWKJ3xiTJ4GBLpKTJfcV82DjxgO0bq3UqbMeuND9oFjBCAhQnH73Ia2paFY3jcPCAnC5gnC5XO4bv8VX8f50xpgCFxBQ8CX+bt3GkZRUkx07nDqkgk38Lk4lfufexKk6/lPKlHE6bUtM9HbPo/5nid8YkydOVU/B3tw9fDjSPZXgPkbBVbmIwKnE78iqxF+2bCkgxOvDShYGlviNMXlSqhReuLmb1jWDc0EpyMTvVPUE41TzOPsvXbr0GeuVKVMaCLISvzHGnM47N3e9nfhL43TPEA10I6vBnUJCnF47rcRvjDGnySrxX3XVAlas2JyPvaY9CdwQqO+lm7ulgRB69eqW5YXFGV3MEr8xxpyhVCkyJf74+Hi+++4qOnd+MB97TbtncBfQtIBL/HCqjv88goOzbpGUn8T/8MNfcfvtX5xdgH5gid8Ykyenl/iT3RX+KSn5aeKZuYTvnTp+gHKEhGSd9tIS/9nU8U+aFMXMmTXPMkLfs8RvjMmT0NAg4uNP3d1NSb8KFM7E7+wqGFgElCU4OOfEf3ZVPUmc6hai8LPEb4zJk+DgwEx99RRM4s+cNL1zczcOKJNtVY8zkDzMmfPtWZT6k7HEb4wptkJCSmXqlvlU4s9POslc4i/Im7uBgYoztm48UJbQ0KwvKmkl/smTX2TQoHvyeBRL/MaYYuzgwX0cOTKCo0ePAhkTfwXOO28yLtfZ9NXvvRJ/UJACoTiJP/c6foBNm3bm8ShJQOkC6VHUFyzxG2PyZPPm9UAFFixYAGRM/PXZvv2+s6wj914dv9OtfyhwEihHaGjuib9Ro6Z5PEoiEExsbNEYwcsviV9EeovIVhH5R0RG+iMGY8zZchJ9Wmue5PTHeJ068pSzerrLeyV+5yHdUCAWp8SfdTWSU8fv3LvIey+dCgRw4sSJs47Tl3ye+EUkEHgLuBJoDPQXkcY5b2WMKTycxH7y5EnnXXqidxK/M5hKXnmvjv9UiT8WCMuljt+5qZuUlNfUGASUYd++ghlDwNv80S3zxcA/qroDQERmAdcCm/wQizEmjypXLs/RozBs2DbGjRvvTvxNgB4ANGv23lmU2C/P9K7gS/ylcRI/hITklPidaqrvv69LjRrj83CUS4DydO48h+rVf8tHtJkFBJxkzpzedO7cucD2Cf5J/LWBPRne78U5a5mIyGBgMMC5557rm8iMMbn64YcXadfuOjp1akFAgFPNs2nTdYSFVWLPnmguuqglp7pg8EyDBmtYtmwsAQHV6dChBkFBHxZYvI880p7AwOdo27Yi69aNpXfvQVmuFxYWxpAhB3j33c506NATEc+bp55zzo+sWvUiTZteRKVKBdeDXUCA5HvA+ayIqnfGzsz2gCI3Ab1U9R73+4HAxao6LLtt2rZtq6tWrfJViMYYUyyIyGpVbXv6fH/c3N0LnJPhfR1gvx/iMMaYEskfif9PoJGI1BeR0sCtwFd+iMMYY0okn9fxq2qKiDwA/AAEAh+o6kZfx2GMMSWVXwZbV9UFwAJ/HNsYY0o6e3LXGGNKGEv8xhhTwljiN8aYEsYSvzHGlDA+f4DrbIjIYeBfoCpwxM/hZMXiyrvCuzRYaQAABfRJREFUGpvFlTcWV975Mra6qhp++swikfjTiMiqrJ5C8zeLK+8Ka2wWV95YXHlXGGKzqh5jjClhLPEbY0wJU9QS/1R/B5ANiyvvCmtsFlfeWFx55/fYilQdvzHGmPwraiV+Y4wx+WSJ3xhjSphCl/hFJFBE1ojIN+73L4nIFhFZLyKfi0hF9/x6IhIvImvdr7d9HNcYEdmX4fhXZVj3SfdA8ltFpJeP45qdIaZdIrLWPd/X52uXiGxwH2uVe15lEVkoItvcfytlWN8n5yybuPz+HcsmrsLyHcsqNr9/z0SkoojMc//bbRaRDoXkO5ZVXH7/jmWiqoXqBTwCfAJ8435/BVDKPT0RmOiergdE+DGuMcCILNZrDKzDGXm6PrAdCPRVXKctewUY7afztQuoetq8F4GR7umRGf4tfXbOsonL79+xbOIqLN+xM2IrDN8z4EPgHvd0aaBiIfmOZRWX379jGV+FqsQvInWAPsB7afNU9UdVTXG/XY4zYpff48rBtcAsVU1U1Z3APzgDzPs0LhER4GbgU28c+yxdi/OfAvfffhnm++ScZaUwfMfyyK/nKyN/fc9EpDzQFXgfQFWTVDUaP3/HsoursH3HClXiByYBj5P9SM13A99leF/fXc3xq4h08UNcD7h/un2Q4SdlVoPJ1/ZxXABdgIOqui3DPF+dLwAFfhSR1SIy2D2vuqpGArj/VnPP9+U5yyqujPz1HcsuLn9/x3KKDfz3PWsAHAamuY/1noiUwf/fseziyshf37F0hSbxi0hf4JCqrs5m+SggBZjpnhUJnKuqrXBXd7ivtr6KawrQEGjpjuWVtE2y2E2Bt5nN7XwB/clcCvPJ+cqgk6q2Bq4E7heRrjms65Nz5pZtXP76juUQl1+/Y7nElsZf37NSQGtgivtYJ3GqdrLjq3OWY1x+/o6lKzSJH+gEXCMiu4BZwOUiMgNARAYBfYEB6q4Yc/9ki3JP/3979/NaxRWGcfz7iIgg3VhE4l67UFHoUqQuLAjdCUWhLhQRFeki/4BIN1J16U7c1La0IkgXgiAIRRBFuij+AEuQUqzULESr4MLEx8WZkJt476gXnBmY5wMhN7kn5OHclzczc27O/EG5ZreuqVy2H9uetf0aOMP8aWNTN5Ovm6+lwE7g17nBDc7X3O97VH2eBi5S5uexpIkq4wQwXQ1vas5G5Wq7xobm6kCNjcwGrdfZQ+Ch7ZvV1xcoDbftGhuVq/UaW6DpRYX3+QC2Mb+IugO4B6xaNGYV1eIM5fTqX2Blg7kmBr4/Sbl+CLCehYtID/iIC2+Lcw3M2e9tzRewAvhk4PH1KtNJFi68nWhyzmpytVpjNblar7FR2TpSZ9eAz6rHx6r6arXGanJ1po/Zbueeux/oNOXFulLWkbhh+xBlAeU7STPALHDI9pMGc52QtJlyuvg3cBDA9l1J5ykv8gxwxPZsg7kAdvP2YluT87UauFi9XkuBn21flnQLOC9pP/AP8DU0Omejck3Rbo2NynWuAzU2NFv1XNt19i3wk6RllEa+j3IVo80aG5XrFh3qY9myISKiZ7p0jT8iIhqQxh8R0TNp/BERPZPGHxHRM2n8ERE9k8YfsYikTwd2S/xPC3fIXDYwTpKu1v2npaRfJK1tJnnE+8nbOSNqSDoGvLB9ashzXwHbbU/W/PwXwB7bBz5eyogPkyP+iPF9A/wGIGmFpEuS/pR0R9Kuasw1YHu1vUFEJ6TxR4xvCzC3Sd4O4JHtTbY3AJcBXPbZmQI2tRMx4m1p/BHjW2n7efX4NuXI/ntJW20/Gxg3DaxpPl7EcGn8EeObkbQEwPZfwOeUPwDHJR0dGLcceNlCvoihct0xYnz3KTsqTklaAzyx/aOkF8DegXHrgLst5IsYKo0/YnyXKFtiTwEbgZOSXgOvgMMAklYDL13dFSqiC/J2zogxVTf6+MH2lzVjJoH/bZ9tLllEvVzjjxhTdRR/5h23ynvK/M2/IzohR/wRET2TI/6IiJ5J44+I6Jk0/oiInknjj4jomTT+iIieeQP43V0mFVUt3AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots()\n",
    "ax.set_title(\"airspeed\")\n",
    "ax.set_xlabel(\"T (s)\")\n",
    "ax.set_ylabel(\"speed (m/s)\")\n",
    "figure1 = ax.plot(ref_time[c:d], indicated_airspeed[c:d],label='indicated_airspeed', c='black')\n",
    "figure2 = ax.plot(ref_time[c:d], true_airspeed[c:d],label='true_airspeed', c='blue',lw=0.7)\n",
    "plt.legend(loc='upper right')\n",
    "\n",
    "plt.savefig('D:/landing_airspeed.png',dpi=1000,bbox_inches = 'tight')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
